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abstract: 

The  determination-,  of  elastic:  properties,  of :  crystals  :by 
sound  velocity  measurements  is:  facilitated  by  use -.of  :'ai 
perturbation  solution,  to  the:  equations  of  motion  for: sound 
waves  propagating  in.  the  crystals:. .  Transformation,  equations 
can  relate  the  solution  for.' arbitrary  porpagat ion  direction 
to  an  established  crystal,  coordinate  system.   A, method  of 
calculating  the  acoustic:  components  of  the  energy; flow  vector 
is  developed  for  several,  high  symmetry  crystal-classes .   The 
validity  of  this-  method  1st  herein-  confirmed  and 'results  are 
presented  which  are  consistent,  with  experimental  1  data.   A 
modern  desk  calculator  was  employed:  to  determine . solutions 
for  hexagonal,,  cubic:  and  tetragonal:  crystals.  . 
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A  major  portion  of  the  thmiL'etUfcHX  development',  presented 
in  this  paper  has  previously,  been:  presented'  by  Prof i. .  John  R. 
Neighbours,  Physics  Department:,.  U...  £..  Nav all  Postgraduate: 
School.   The  mast  recent  paper,.  co:-authored.  by  Assoc iate^ 
Professor  G.  E.  Schacher ,.  Physics:  Department-,,  Ul. .  S.~, .  Naval  1 
Postgraduate  School,  is  titled  "Determination:  of:' Elastic . 
Constants  from  Sound  "Velocity  Measurements"  in:  Crystals:  of  : 
General  Symmetry."'   In  addition:,,  several,  investigators:  have 
published  methods  similar  to  this:  development-.. .  These:  papers 
are  included  in  the  List  of  References:.. 


I.   INTRODUCTION 

The  task  of  relating  velocity  measurements:  tax  elastic: 
constants  of  a  crystal  is  complicated  by  tire  fact  that",  the 
eigenvectors  which  describe  the  direction,  of  vibration  in 
response  to  the  sound  wave  are  neither  perpendicular"  nor 
parallel  to  the  propagation  direction...  The  energy  flow  may, . 
however,  be  described  by  modes  of  vibration,  which  are:  quasi--- 
longitudinal  (essentially  parallel  to  the  propagation. direc- 
tion).  In  practice  one  vibrational  mode  dominates  the  others: 
for  a  given  propagation  direction..   In.  addition,  the:  transducer: 
used  to  measure  the  velocity  has  greater  response  to:  one:  mode: 
of  vibration  so  it  is  possible  to  separate  modes  by  appropriate: 
choice  of  transducers. 

For  certain  high  symmetry  crystal  classes,  special!  direc- 
tions of  propagation  exist  in  which  the  response  "is  pure 
longitudinal  and  pure  transverse ..  Velocity  measurements  along; 
these  directions  can  be  directly  related  to  crystal,  elastic 
constants.   However  the  energy  flow  is  not  necessarily  parallel- 
to  the  propagation. 

For  either  high  or  low  symmetry  cases  the  elastic  constants: 
can  be  determined  by  solution  of  the  crystal  vibration  equation 
of  motion.   This  process  may  be  simplified  by  appropriate 
choice  of  coordinate  system  and  application  of  several  assump- 
tions.  The  resultant  secular  equation  is  solved  by  a 
perturbation  method  and  the  values  for  the  elastic:  constants: 
may  be  refined  by  iteration.   The  constants  thus  obtained  are 


linear  combinations:  of  the:  elastic:  constants: defined :by ; 
conventional   coordinate   systems:  and:  can.bee converted  :by 
rotational  trans^oizraatijan-..     The:  angle-:  of:' deviation: of  :"t he 
energy  flow  vector  from  the  propagation  direction   is : of : 
particular   interest:  in.  order   to.  construct:  a:-,  physical  Ipicture 
of  the   crystal  response-..     This:  angle:  can  be:  obtained  : by 
application  of  the  general  energy   flow  equations :  to  : the. -. 
coordinate   system  defined,  by  the.  propagation,  direction.  . 


II.   SECULAR  EQUATION 


The  general  form  of  Hooke's  Law  is 


ij    ijkl  kl 


(1) 


Two  common  conventions  of  subscript  notation  for  the  above 
equation  are  shown  below. 


^r 


FIGURE  1 


The  strain,  S,  is  related  to  displacement  by  the  equation 


Skl   =  2  (  Uk,l  +  Ul,k) 


.    3u,    3un 

-  (  — -  +  — -  )    (2) 
2  k  3U,    9uk  ;.    K    J 


where  u  is  displacment.   In  contracted  notation  Hooke's  Law 
may  be  written  as 


T.  =  C. .S, 
1    ij  j 


where  the  following  notation  applied: 


(3) 


X   =  T..,  =  T, 
x     11     1 


Ty   T22   T2 
Zz  =  T33  =  T3 
Yz  m  S  '   r23  =  T32  *  =*  W 

Xz  "  Zx  "  TL3  =  T3I  "  % 
y     x     12     21     a 

The  existence  of  a  strain  energy  function,  reduces:  the:  3.6. 
elastic  constants  to  21  independent  constants:  for.  the:  most 
general  case.   This  number  can.  be  further  reduced,  for  specific 
crystal  calsses  because  of  symmetry  considerations:. . 

If  an  arbitrary  coordinate  system  jis  chosen,  such,  that  the 
x'  axis  lies  along  the  direction,  of  propagation  the  equations 
of  motion  can  be  written  in  terms  of  the  elastic  constants  in 
the  arbitrary  coordinate  system  and  displacement:..  In:. this: 
development  uT,  vr  ,  wr  are  displacements  along;  the:  x' ,.  y;' ,  . 
and  z'  axes  respectively-  The  primed  notation  emphasizes  the 
fact  that  choice  of  coordinate  system  is.  arbitrary.   The 
equations  of  motion  are 


2  3T     3T     3T 

Fx,  =  ma  =  p  2-gi  dx'dy'dz'  =  p  (—£  +  --$  H-  — ?  )  dx'dy'dz' 


2 

3  U     3 

^2~  =  3^  (  C11S1  +  C12S2  +  C13S3  +  Cl4Si  *  CV5S5   *  Cl6S6 

iTc         3T,. 
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If  plane  wave  pirapagatijon.  is:  assumed:'  ther  partial  1  deriva- 
tives with  respect  to  y'-  and:',  z:''  vanish,  andi  thef  equation 
becomes 

2     .  Z  2" 

11   a*-2  *  3X-2'  «-=  a*-2:  (6) 

Similarly,   the  equations  of"  mot±ori- TiDr/ the:  remaining  ^directions 
are 


F  -  a   *V      =.     CT"      ^'  Id     lB±l+-C'      i^- 

*'  P    *t2  ^3x.'Z         ^3^2  ^:ax,2 


An  additional  assumption,  is:  that",  the:  solutions:  are  .periodic 
of  the  form 

uf'  -  A -j-  cos  (kx'-wtr.) 

\rr  =  Az   cos  (kx'-wtr)  (7) 

wr  =  Ao  cos  (kx'— ujt) 

The  assumed  periodic  solutions  are  inserted  in  the  equations 
of  motion  and  the  partial  derivatives  are  evaluated.   The 
resultant  equations  may  be  further  simplified  by  the  relation 
V  =  w/k  where  V  represents  velocity  and  is  not  to  be  confused 
with  v',  the  displacement  in  the  y1  direction.   The  equations 
may  now  be  written 


il  -  rl 

ftZClfi 

A    C  ' 

Vie 

%Cc&  -  vZ) 

A3Cf56 

Aici5 

*ZC& 

%(C55  ~V 

21 


u 


V* 


w 


=   0  (8) 


The  secular  equation,  has:  three:  rea-X.  roots  corresponding  to 
one  quasi-longitudinal  and  twa  transverse  modes: of : vibration. 
If  the  solutions  are  assumed:'  toe  he  essentially  longitudinal 
or  transverse  with  small  deviations  from  these-:  directions 
the  secular  equation,  yield's:  a-,  perturbation  solution:: 


°11 


=  P 


V?  - 


r 


a, 


2 

a  a  r ' 


ALA2Cj| 


|[A3(C11-Cfe):   A2(Cil  -C66) 


Note  that  the  amplitudes  of  the  assumed  periodic:  solutions  do 
not  affect  the  caXculation..  The:  remaining  solutions  to  the 
perturbation  equations  are. 


C66 


=  *s- 


c2 
Cl6 


'56 


[L66   uIi 


°6.6l'   °55 


C55 


=   pV|  - 


c2 
°15 


,,2 

'56 


L55    11 


L55    66 


U 

i 


(9) 


In  particular,  if  the  terms  in  square  brackets  vanish  the 
solutions  correspond  to  pure  mode  response  and  the  elastic 
constants  are  determined  directly  from  the  measured  velocities. 
When  the  square  bracketed  terms  are  present  velocity  measure- 
ments along  several  directions"  yield  an  adequate  number  of 
equations  to  determine  all  six  elastic  constants  of  the  secular 
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equation  to  first  order..  Tftessase  values:  aasee  then,  inserted  in 
the  square  bracketed  terms  and:  the:  solution r  is:  refined.   The 
primed  elastic  constants  are:-  linear  combinations:  of:' those 
defined  by  the  conventional-  coordinate:-  system.   The  rotational 
transformation  relating  the  two  systems:  is;  developed  in 
Section  IV. 

The  amplitudes  of  the  periodic.  solutions r may- be ^determined 
from  the  equations:  cf7  motion...  I1V  the:  amplitude:  of :" one  mode 
of  interest  is  fixed  at  an.  arhitrary value  the: remaining  two 
amplitudes  can  be  expressed  in  terms:  of  _"  the  elastic  .constants 
and  the  measured  velocities:..  For.  example  ,.  if :  the  amplitude 
A,  is  set  equal  to  unity  the:  longitudinal^  mode" ■-  solution  is 


g4  -  4i 

Ci5 


4.5= 


'16 


45 


L      46 


^  '  C55 


C56 


oV3:--c55 


'4  -  4i 

46 

+ 

°3s 

C56 

cV3   -  C55 

+ 

45- 

C16 


(10) 
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ENEFKIY"  FLOW 


The  components  of  the:  Energy  Flow  vector,  are- ' 


E  .  =  urXr  +  v'X'  +  w'X1 
xr       x      y       z: 

E  .  =  urYr    +■  vT1"  +■  w'T1  (11) 

yr      x      y      z: 

E  .  =  urZr  +  v^Z"'"  +-   wrZ~' 
zr      x      y      zl 


The  x*  component  of  this  vector  can  be  calculated:  as: follows 
Ex»   =  u"  C  C1I  aF""  +"  C1E   a^  +  C15  3*T* 

*  v  *  C1.6  ax'  *  C6.6.  3x'  +  °56  3x:'  ;      U^; 
4  v  (  C15  9x'  +  %'  3x'  +  C55f3x'  J 


Note  that  the  plane  wave  assumption  has  been  applied -to  yield 
only  those  terms  which  contain,  partial  derivatives: with  respect 
to  x'.   The  transverse  components  of  the  Energy  Flow  vector  are 

+  ;'  <C12  IF"  +  C26  Wr   +  C25  IP-  )       (13) 
+  w  cci4  ax'  +  c46  ax'   c45  ax^  } 


A.  E.  H.  Love,  A  Treatise  on  the  Mathematical  Theory  of 
Elasticity  (New  York,  IpTJT  ch.  7. 
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E  .   =   u»  (  C'   |Hl  -f  C^      |Zl  +  ar   |^1  ) 

z1  15  3xr     66  ax'     5:5  3-x:1 


+  v'fC   iHl  +  rr   iZl  +  Ct-  Mil) 


'14  3xr    4-6"  5x:r    Jl-5"  3x.' 


+  w   C  C13  3x'  +  L36  3x'  *  G25  3x'  ; 


On  substitution  of  the  periodic  solutions:  the:  x:' '  component. 
becomes 

Ex'   =   wAisin  n  [Ci1(-kA1sin  n )  +  C-T^C-kA^siir.  n ) 

+  C'  (-kA^sin  n )  ]   +  ajA^ain:  n  CG^'-kA^.  sih  n) 


+  C^6(-kA2sin  n )  +  C^C-kA  sin  n)I 


+  uA^sin  n  [ ] 


where  n  =  kxf  -  wt . 


The  above  form  may  be  further  simplified  by  defining. 


x'  2 

-wksin  r\ 


The  final  form  of  the  Energy  Flov;  components  is  shown  below. 
Px,  =  C^A^  +  c<66k\   ♦  C'5A2  +  2Ci6AlA2  +  ZC^AjA   +  2C^*2A3 


P 

y 


•  =  Cl6Al  +  C26A2  +  Ci5A3  +  (C66  +  Ci2)AlA2  +  (C56  ^CVA1A: 


(C-5  ♦  Ci6)A2A3 
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V  =  Ci5Al  +  C46A2  +  C35A3  +  (C56  +  Cl4)AlA2  +  (C55  +  Ci3)AlA3 


(C45  +  C36)A2A3 


(15) 


The  resultant  Energy  Flow  vector  determines  deviation  from 
direction  of  propagation  as  shown  in  the  diagram  below.   Note 
that  the  angle  5  is  the  compliment  of  the  polarization  angle. 


Px' 


z»' 

**«*    P 

7,' 

^^<^      ,py'    y1 

5^  i           1     /'   * 

k 


FIGURE  2 
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IV.   TRANSFORMATION  OF  ELASTIC  CONSTANTS 

The  desired  elastic  constants  are  defined  by  a  fixed, 
unprimed  coordinate  system  which  is  related  to  the  crystal 

basis  vectors  a,  b,  c  in  accordance  with  convention  adopted 

2 

by  the  I.R.E.  standards  committee.   Specific  examples  of  this 

convention  are  provided  in  Sections  V,  VI  and  VII.   The  primed 
coordinate  system  is  chosen  so  that  the  x'  axis  lies  in  the 
x-y  plane.   This  leads  to  a  relationship  between  the  two 
systems  in  terms  of  the  direction  cosines 


£  =  sin  e  cos  <\> 
m  =  sin  6  sin  <j> 
n  =  cos  e 


(16) 


FIGURE  3 


Standard  Committee,  "Standards  on  Piezoelectric  Crystals," 
Proc.  I.R.E.  37  (September,  19^9),  p.  1378. 
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Tbe  transformation,  elements,  are   given:'. Inr. tabular:  form 
below. 


X. 

y 

z: 

X1 

£ 

m 

nr. 

yr 

-1 
-ma 

£a 

o: 

ZLV 

'.  — inct~ 

i 

-1 

-mna 

a: 

(17) 


The  elements  nave  been  reduced:  by  use  of:  the:  definition 

2  _   2    2 
a  =    I      +  m  . 


For  example: 


x  —  £x:r   —ma   y'  -Una.     z:1  ' 


(18) 


The  transformation  of:  elastic  const  ant  s.may"  be  :  obtained 
from  Hookers  Law; 


%  "  °«*j 


(19) 


If  the  trans  format  ions  of/  stress  and  strain:  are  -TT  =  6T. 
and   Sr  =  yS~      the  elastic:  constants  are  related  :.by  the 
expression 


Cij   -  6(V 


-1 


(20) 


where  y 
.2 


2  -2 

m  a 


-1 


=   6   and   8  is~  the  matrix  shown  below. 


m 

n2  -2 
l    a 


2  2-2    2  2-2 

I    n  a     mna 


-2 

2.  mna 

2   -1 

-i   na 


-2 
-£mna 

2   -1 
-m  na 


n 
0 
a 
0 


2mn 
0 
-2mn 


2£n 
0 
-2  Jin 
-nr 


2£m 


-2£ma 


-2 


„   2  -2 

£mn  a 

/02   2.  -2 
-n(£  -m  )a 


-£ma 


-1 


£ma 


-1 


na  m(l-2n2)a 
-1 


£(l-2n2)   1_  -2£mna  -1 


-mna 


-1"_     ,  „2   2s  -1 
(£  -m  )a 


0       £na 

The  Energy  Flow  equations  and  the  rotational . transformation 
may  now  be  applied  to  specific  crystal  classes. 
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V.   CRYSTALS  OF  HEXAGONAL  SYMMETRY: 

In  application  to  certain,  high:  symmetry  crystal! classes  : 
the  method  Is  simplified  due  tn  symmetry  considerations :and : 
the  fact  that  certain  primed  elastic:  constant's:  may  vanish. 
This  leads  to  a  partially  diaganaX.  secular.' equation:  and: the 
existence  of  semi-pure  modes-  of"  vibration. .  For:  hexagonal  1 
crystals  the  I.R...E..  convention  defines:  the:  unprimed:  coordinate 
system  such  that  the  z-axls  Is:  para 11  ell  t"o:  the:  crystal ^basis  : 
vector  c,  i.e.  the  sixfold  symmetry  axis:, .  and:  the:  x-axis  :  is  : 
parallel  to  a,  the  two  fold  rotation:  axis:. .  The:  y-axis  :is  : 
perpendicular  to  x  and  zr  to  preserve:  a  right:  handed: coordinate 
system.   With  this  definition  the.  primed:  elastic:  constants  . 
become 

Cll  =  a^Gll  +  2nZa?(cx3-  +  2n:l^)    "*"  n  c33 

C£5  =  n^a2^  -  nZaZ(2C:LI  -  C^)  +- (i:-2h.2")2C44 

He  =  lh2(cii  ~aizy  +  n\*  (22) 

C'   =  -na3Ci:L  +  na(I-2n2)  (C"^.  +■  20^)  +-  n^aC^^ 


C16  "  c56  "  ° 


Since  C',-  and  C'^  are  zero  the  existence  of  a  semi-pure  mode 

Is  Implied.   Physically  this  means  thatr  the.  vibrational "response 

of  the  crystal  lies  along  one  of  the  primed  system  axes  for. 
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one  mode  (pure  transverse)  and  the  remaining  modes  are  quasi- 
longitudinal  and  quasi-transverse.   The  relative  orientation 
of  the  eigenvector  components  to  the  propagation  is  shown 
below. 


3  -  (  QT) 


V2  -  Pure  transverse 


y' 


FIGURE  4 

Note  that  if  C',-  were  to  vanish  the  secular  equation 
would  become  diagonal  and  the  response  would  be  pure  longi- 
tudinal or  transverse.   This  situation  may  occur  for  materials 
possessing  hexagonal  symmetry.   The  existence  of  this  acciden- 
tal  pure  mode  was  predicted  by  BorgnisJ  for  the  case  of 
propagation  in  a  direction  that  makes  an  angle  8  with  the  c 
axis  defined  by 


2      °33  "  °13  "  2C44 
tan   0  -  p     p ^p 

ull  "  13  "   44 


(23) 


The  above  requirement  corresponds  to  the  fact  that  C'   vanishes 
and  the  secular  equation  is  diagonal  in  the  perturbation 
solution. 

Equations  (22)  are  sufficient  to  calculate  the  velocities 
of  each  of  the  modes  from  known  values  of  conventional  elastic 


F.  E.  Borgnis,  "Specific  Directions  of  Longitudinal  Wave 
Propagation  in  Anisotropic  Media,"  Phys.  Rev.  98  (August, 
1955),  p.  1000. 
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constants.   The  elastic  constants:  cdT  numerous:  hexagonal! 
crystals  have  already  been,  determineil  and:  the:  application-  of. 
the  perturbation  method  may  be  demon strated:.  by  comparison,  of  ' 
calculated  mode  velocities  with  measured:  va-lues:. .  The:  solu- 
tion for  energy  flow,  however,,  requires  nine,  additional- 
constants  i.e.  C-^2r  ^13''  ^T4-''  ^25:3:     eCS>:   ^355'  ^36° ■'  ^45' :  ^46  ' 
Some  of  these  may  vanish,  as  a  resuit  of:  symmetry..  For: 

hexagonal  symmetry  C''^,,.  CAV,.  C^  and  C"''^.  are  alll  identically 
zero.   The  remaining  primed  constants  are- 

Ci2    =   A12    *  AI3 

c«     =  Ci  -  2n2a2)cL2  +-  nza?(:a:ix  +-  c33-  —  4C^_)' 

C^    =   na(C13    -   C12)  (24) 

C'      =   -naCn2^  +    ( l-2"nZ)  C^  -  0^33,  +*  2i(i:-2h2-)Cij1| ) 


C46  =  "  |n^Cn  "Ci2  -bzw) 


In  addition,  the  fact  that  C'V  and"  CT'.,-  are  zero:  leads:  to 
simplified  forms  of  the  amplitude  equations.   With  these 
simplifications  the  equations  for  energy  flow  in  hexagonal 
crystals  are 

Px'  ■  CilAl  +  C66A2  +  C55A3  +  2Ci.5AlA3 


V    =    (C66   +  Ci2)AlA2   +   (C25   +  Cie)A2A3  (25) 


Pz,  -  CJgif  +  C'^k\   *    {C<55  *   C-3)A1A3.  t-  C^A2; 
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Some  of  the  terms  in.  the-  ahove-  equations:  will!  vanish.- when 
certain  roots  of  the  secular-  equation',  are:-  applied  .  since-  the 
amplitudes  of  the  assumed,  periodic  response:  vanish.. ,  In.  the 
notation  used  below  the  second  subscript,  denotes'  the? root" of 
the  secular  equation  that  is  used  in  the:  calculation. .  For 
instance,  A,,  is  the  amplitude  when,  the"  root.  V-.  -  is.  used. 
The  amplitudes  are 

pa£— C*' 

A    =  7  &    —  CT    A:    —  — =^~ 

All  ~  •  ^12   u   R13-      C'  - 

l£>; 


A21  =  Q  A22   =  L   A~:  =  0:  (26) 


t 


Pvf  -  q 

a    -  ___± ±i.    /v    —  0"    A- 

A31     c^     *1Z      u   A33.: 


The  energy  flow  relationships7  for  the-,  modes,  of :" vibrational 
response  in  hexagonal  materials  are 

For  the  root  V,-  (Quasi-LongitudinalO 

V  "  c£lAII  +  C55A3X  +  SCI5A11A31: 

P   ,  =  0  (27) 

y 

V  "  CL5AL1  +  C35A31  +  (C55  +  CI3)A11A31. 

For  the  root  Vp  (Pure  Transverse) 

2 

V  *  C66A22 

Py,  =  0  (28) 

2 

V  "  C4oA22 
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Far-  the  root  T~  ( Qu  a:ST-Tr_ans.y  erse^) 

2-  7 

*xr   uIiA13  *  u55  331     15:  13- 33; 

Fyf  -  Q  (29) 

Pz''  =:  ^I^H*  C2^t32  *"  (C55f  "*"  C13:)A13A33: 

For  tills  symmetry  the  above  three  equations:  (27",  28 ^  .and  29) 
show  that  the  yr  component  of  the  energy  flow  vector. is 
always  zero.   Thus  the  angle  \p.   defined"  on.  page  111 -is:  always 
zero  and  any  deviation.  of  the  energy  flow  vector  from  the 
propagation  vector  takes  place  in  the  x:'--Z-  '  plane i  . 

In  order  to  demonstrate  the:- validity- of:"  the:  perturbation 
method  solutions  were  computed.:  on  ai  Hewlett  Packard:'  9810A 
desk  calculator  and  compared  to  published:,  data:  for:  hexagonal 
crystals,  specifically  zinc:  and  beryl!,,  as: presented. by  M.  J. 
P.  Musgrave..   Transformation  equations:  using  known/values  of 
elastic  constants  were  programmed  and  the:-  velocities  obtained 
were  in  agreement  to  ..G5X.-  The:  energy  flow  components  were 
calculated  and  the  angles  of"  deviation  from  propagation 
direction  were  consistent  to  within  2°    of  those  presented  by 
Musgrave.   The  graphs  on  the  following  pages  are  examples  of 
results  obtained  for  x-z:  plane  propagation  in  hexagonal 
crystals.   In  these  graphs  the  deviation  angle,  6,  is  plotted 
as  a  function  of  the  angle,.  6,.  which  the  propagation  vector 
makes  with  the  z-axis  for  each  of  the  three  modes.   In  addition 


4 

M.  J.  P..  Musgrave,  "On  the  Propagation  of  Elastic  Waves 

in  Aeolotropic  Media,"  Proc .  Roy.  Soc .  226  (195*0  p.  356. 
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the  sign  of  all  deviation,  angles  in.  the:  cubic:  and:, hexagonal 
cases  has  been  reversed  fir.  order"  to:  compare:  the:  results  with 
Musgrave's  data.   Table  1  lists,  the:  known,  elastic:  constants 
and  densities  required  far  the  aaicul  ation. 


TABLE    I 
Elastic  Constants 

('10        dynes/era    ) 

Density 
(gm/crrr) 

I    cn 

aiz 

C13- 

C33: 

ci|i|. 

Zinc 

Magnesium 

j        14.3 

6*..  35 

2E..59F 

3:.  3: 
2.  XT 

5:.o: 
6:64- 

i:8«- 

7.1 

1.75 

The  metals  Zinc  and  Magnesium  are  two  different: examples 
of  hexagonal  material...  The  difference  between. C,,- and: C-~ 
in  the  two  materials  affects  the  response:  characteristics . 
In  both  materials  the  Pure  Transverse  mode  (TO  has ;.  a:  non-zero 
deviation  angle  for  energy  flow..  This,  means:  that:  although 
the  mode  is  a  pure  vibration,,  associated,  with  the: root  Vp  and 
lying  along  the  v'  axis,,  the  energy  flow  is,  non-zero  and  not 
coincident  with  the  x'  axis  (the  propagation  direction). 

The  response  for  Zinc  indicates  that  no  accidental. mode 
exists  since  the  Quasi-Longitudinal  (QL)  response  does  not 
become  zero  except  at  the  end  points  8  —0°  (the  c  direction) 
and  8  =  90°  (the  x-y  base  plane).   The  absence  of  accidental 


W.  P.  Mason.,,  ed .,.  Physical  Acoustics  v.  HIP/  (New  York, 
1965),  p.  32. 
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response  can  be  verified  by  equation  (23).   The:  Quasi-- 
Transverse  (QT)  mode  has  a  zero  value  between. the vend .points 
but  this  does  not  imply  pure  made:  response  since:  the ;  amplitude 
A,-,  has  a  non-zero  value  at  this  point.   For.  Zinc:  the- general 
effect  of  elastic  anisotrapy  on  the  energy  flow  is ~ to  channel 
the  energy  of  most  of  the  waves  toward  the  base: plane. 

An  accidental  pure  mode  Is  observable  in  Magnesium.   The 
Quasi-Longitudinal  mode  has  a  zero:  value  at  8  ==51*6°. 
Analysis  of  the  amplitude  shows  that  both  A,-,- and:. A~, .  are 
zero  along  this  direction,  of/  propagation.   Thus:  the  amplitudes, 
and  hence  the  vibrational  responses  are  colinear:  with  the 
primed  coordinate  system..  The  Quasi-Transverse:  mode  :  has 
essentially  the  same  form  of  response  as  Zinc: but  is. opposite- 
ly directed.   This  Is  due  to  tha  difference  in. the  -  value  of 
(C,,  -  C_^)  for  the  two  materials..  The  magnitude: of  the 
deviation  angles  In  Magnesium  Indicates  that  the  material  is 
essentially  isotropic 
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VI..   CUBIC:  SYMMETRY' CRYSTALS' 

Calculations  for/  cubic:  symmetry-  crystals:  follow  the  same 
pattern  as  that  for  hexagonal,  crystals:. .  The:  conventional 
axes  are  parallel  to  the:  crystal!  basis:  vector s: a, .b ,  and  c. 
The  elastic  constants  are:  (X,---  "  C'-,, .  C,  -,  -  ==  C,~,  .and:  C/-g  =  C^h. 
Transformation  of  these  constants:  leads:  to:  the -equations 
shown  below.. 

h  II  h.  P"  ?"    ?'  ?'    p  '  ?  ' 

C-j^  =  (2   +  m  4*- rr.  yCL-j-  +-  2U  m  +- ■£  fi  +-m n  .)(C12  +  2C4ij) 

C'   «  2n2a~2U   +- m  +  iSnr^XC^-  —  C12-  —  20^  )  +  C^ 


(30) 


2.  2  —2 

clfi  =  2*  m  a     ^cn  ""  ci2:  ~~  2;c4ip  +'  ci\i\- 

C<      =  -nxx~L(£   +  ra  -  n^MC^  —  C12- —  2C^) 


For  propagation  in  the  x— y,.  x-z:, .  y--z:  or.  (110)  planes  the 
constants  CJV  and  C'g  vanish. yielding  a  partially  diagonal 
secular  equation  as  for  the  hexagonal-  case.   The  velocity 
calculations  are  identical  to  those  for  hexagonal  materials. 
The  relations  from  the  secular  equations  are 

V2   =  C66/P 


Vl,3  =  27  ((Cil  +  C55}  ±  ^/('Cil+C55)2  -4<cil?55  "^    } 

(3D 

Further  response  analysis-  in  the  above  planes  yields 

amplitudes  that  vanish  for  certain  modes  of  propagation  and 
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elastic  constants  that  vanisli.  toy  choicer  of.  pd'ane.  of  ^propaga- 
tion.  Thus  only  Chg-*  C.L  and:  CL1'  remain,  to:  dbt  inserted  :  in 
the  resultant  energy  flow  equations:  which  are-,  similar:  to 
Eq.  (25)  for  hexagonal  symmetry..  The.  graphs:  on  the:  following 
pages  show  solutions  for  Aluminum  with  sound:  velocity;  propa- 
gation In  the  (1QQ)  plane  (x-y  plane.)  and  the:  (11D)  plane 
((j,  =  45°).   The  results  agree-  with  Musgrave:' s:  data.  -  The- 
elastlc  constants  far  Aluminum  are 

11  2' 

C,,  -  5. .5:   x  Id    dynes/cm 

XX  2 

C-j-o  =  H--9F  x:  10    dynes/cm 

XT  2' 

GTjih  —  Z..o   x.  10    dynes/cm 

The  value  for  density  is  2L..6.5"  gm/cm  . .  These:  values: are : taken 
from  Musgrave . 

The  results  for  (100)  plane  propagation  in. Aluminum- show 
the  existence  of  the  accidental  mode  at'  6  =-  45D:  which  is 
predicted  by  Borgnis..   The  high  symmetry  off  ther  response-. 
curves  is  expected  In  the  case  of  cubic:  crystals- .  In: addition 
the  deviation  angle  for  the  Pure  Transverse  mode-  is : zero. 
Thus  the  energy  flow  Is  colinear  with  the  propagation  direc- 
tion.  For  propagation  in  the  (110)  and  (100)  planes  the 
angle,  i> ,  is  zero. 

The  response  for  (110)  propagation  merits  particular 
attention.   The  Quasi-Longitudinal  mode  deviation  angle  goes 
to  zero  at  6  =  5^.6°  which  is  the  [Hi],  direction.   This  is 
not  an  accidental  mode  but  is  Indicative  of:"  the:  fact": that 
C!,-  is  zero  because  of  propagation  direction.   Thus  the 
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secular  equation,  is:  diagonal,  and.  the  three: modes: are  pure. 
In  addition,  the  velocities-  Vp-  and  V-  are. equal land : the  devi- 
ation  angles  for  the.se.  manfess  are:  equal  and.  non-zero.   The 
phenomenon  of"  internal-  cccnicai  refraction  exists: for  this 

propagation  direction..  It  can.be  shown  that  P:  ,  . and 

2     2   1/2 

Cp  r  +  P  „ )    are  constant. .  Since  the  response  for  A-,  is 

equal  and  opposite  to  that,  for  A?,  the  response  ;  for  any- 
linear  combination  of  the-  two:,,  (making  an  angle  y,  .with  the 
yr-axis)  would  occur  at.  twice:  the  angle  Y>.in:the:-y  direction 
In  other  words  the  energy  flow  describes  a  cone : about  the 
propagation,  direction..  If"  the:  amplitude  fbr:  one  mode  is 
rotated  the  response-  rotates.-  in  the  opposite  direction  at 
twice  the  angular  frequency  of-"  rotation. 
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VH.   TETRAGONAL'.  SYMMETRY. 

Solutions  for  tetragonal,  crystals:  represent":  a-  step: 
forward  since  there  is  little-  data,  avai  Table::  to:  be-:  compared 
with.   The  I.R.E..  conventions"  define:  the:  z:— axis:  parallel! to 
the  basis  vector  c  and  the  x— axis:  paraTlJell  to:  the--  basis: 
vector  a,  or  the  twofold  symmetry  axis-..  The-- two:  types  .  of : 
tetragonal  class  materials  are  distinguishable-  by-  the^ fact: 
that  C-,g  =  0  for  Tetragonal  E  materials:. .  The:- transformation. 
of  elastic  constants  is  otherwise  identical",.  The:  transf orma-- 
tions  for  Tetragonal  II  materials"  are: 

Cll   =    ^   +  m^cii  +  Zs^vfajz  +  ar^aSc^  -H  2E^)    +-n?C~    , 

+   4z2m2Cg6   +  4£mU2  -  ra^C^ 

C'      =  n2a"2U4   +  m^)cn  +   28^Frr?cc^C^2  *  2C66) 

-  n2a2(2Cl3    -   C      )    4f  (X  -  ?nZ)ZC^  4-  4-  mn"2a~2{  £2   --m2  )Cl6 

C^6    =    2£2m2a~2(Ci;L  -  CL^)    +-  nZC^  +- -  aT2^?  —  nr2/)2C66. 

-  4£ma"2(£2   -  m2)Cl6  (33) 

C'      =   -ncT1^    +  m4)Ci:L   -   2£2m2n   "^C^-  +  naill  —  2n:2MC£3    +"20^) 

o  2    2      —1  1      2  2 

+  n  oCqq   -  4fc  m  na     Cgg   -   4ij,mna      (£      —  ra   )C£g 

C^_6    =    -£ma_1(5,2    -   m2)(Ci;L   -   C12   -    2C66) 
+   a    ■L(£^    +  m      -    6im    )Cl6 

C'6    =    £mna"2(£2    -  m2)(Ci;L  -  C.I2   -  ^55) 

-2,    4  4'        Cfl.Z  22Vn. 

-na       (£      +  m     —  oJl   m    )C,.g 
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Once   again  the   ctLai.ce  aT   (XQOL)    car  (ITO.)   planee  propagation 
eliminates   the  constants-  CjL*  and   (Xj-V  and:  the:  secular:"  equation 
is   partially  diagonal..      The  energy   fiXaw  calculations:  are :- 
simplified    since   GXh*   Ca6>-  c\.q  sn^L  ^/ict  ar-er  identically-  zero. 
The   remaining  constants  are 

C'12   =    (l-n2)Cl2   +  ei2C12 


(33) 


C^   =  n2(l-n2)(C11   +  C22  -4C44)    H-  (X-2h^(^Hfi21)C^ 

C^   =   -n(l-n2)1/2    (CL3  -  CI2) 

C'      =   -n3(l-n2)1/2   CTT   -  n:(i-n2:)I/2(l-Zn2:)Cn^ 
3p  il  lo: 

+  n(l-n2)3/2C33   -  2n.(T-nZ)X/2L(±--2n2:)CnA 
Ci6   =  n(l-n2)1/2(C^   -   C^) 

The  resultant   velocity  and  energy-  daw  calculations: are: 
identical   to   those   for   hexagonaX  symmetry..     The:  results: are 
shown   on   the   next   graph  for    (I0TO-)    plane"  propagation  .in:  Tin. 
The  required  data  for  calculation   Is   given. below. . 

C1X   =   8.27   x   1Q11  dynes/cm2" 

C12  =  5.78 

C13  =3.^2  " 

C33   =10.31 


C41|   =   2.69 
C66   =   2.82 


=   7.39   gm/cm- 


6Ibid. 
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The  response  is  gxmt'T &r   tec  that.  £exe  Zinc:  iir. that"  no: 
accidental  mode  occurs..   In.  ad'dditidn.  the:  small! deviation 
angles  again  indicate  that  the:-  material!  is:  essentially; 
isotropic.   For-  the  tetraganai  results.  the_  deviation, angles 
are  not  reversed  in  sign. 
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vni..  conclusion 

Agreement  with  published:  data:-  indicates:  that:  the  method 
presented  In  this  paper  is  a  valid:  one:  and:  is:  particularly 
useful  in  the  determination:  of  energy  flow  components .   Thus 
it  is  possible  to  calculate  the  direction,  of:"  energy  flow  for 
waves  propagating  in  an  arbitrary  direction. .  The: results' 
also  show  that,  perhaps  surprisingly,,  there:  may-  be.  directions 
of  propagation  for  which  all  the:  modes:  are:  pure  but"  the  trans- 
verse wave  energy  flow  vector  and:  the:  propagation  vector  are 
not  colinear.   This  method  oan.be  extended:  to:  a  crystal: of 
general  symmetry  ►   The  effects  off  crystal,  stability  and  .multi- 
ple mode  response  are  areas  for  further  investigation. 
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13.  ABSTRAC  T 


The  determination  of"  elastic:  properties  of  crystals  by  sound 
velocity  measurements  is  facilitated-  by  use  of  a  perturbation 
solution  to  the  equations  of  motion  for  sound  waves  propagating  in 
the  crystals.   Transformation  equations"  can  relate  the  solution  for 
arbitrary  propagation  direction  to:  an  established  crystal  coordinate 
system.   A  method  of  calculating  the  acoustic  components  of  the 
energy  flow  vector  is  developed  for  several  high  symmetry  crystal 
classes.   The  validity  of  this  method  is  herein  confirmed  and 
results  are  presented  which  are  consistent  with  experimental  data. 
A  modern  desk  calculator  was  employed  to  determine  solutions  for 
hexagonal,  cubic  and  tetragonai  crystals . 
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